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DEFINITION  OF  SYMBOLS 


p  *  Total  pressure  (mb) 

E^  *  Saturation  vapor  pressure  of  moist  air* 

Ew  «  Saturation  vapor  pressure  (over  water) 

Ej  *  Saturation  vapor  pressure  (over  Ice) 

E$  *  Saturation  vapor  pressure  (reference  unspecified) 

Er  *  Saturation  vapor  pressure  -  Goff  Gratch  formulae 

“  (reference  unspecified) 

Et  *  Saturation  vapor  pressure  *  TetSns  formulae 

(reference  unspecified) 

ET1  *  Saturation  vapor  pressure  -  Table  look-up  procedure 

IL  (reference  unspecified) 

E  »  Saturation  vapor  pressure  -  Polynomial  formulae 
L  (reference  unspecified) 

Dr  =  Derivative  of  saturation  vapor  pressure  with  respect 
"  to  temperature  -  Goff-Gratch  formulae  (reference 

unspeci fied) 

Dt  *  Derivative  of  saturation  vapor  pressure  with  respect 
to  temperature  -  Tetens'  formulae  (reference 
unspeci fied) 

D.  =  Derivative  of  saturation  vapor  pressure  with  respect 
L  to  temperature  -  polynomial  formulae  (reference 

unspecified) 

c  *  Ratio  of  molecular  weights  (.62197) 

r^  *  Mixing  ratio  of  water  vapor  (water  reference) 

t  *  Temperature  in  degrees  centigrade  (Celsius) 

T  *  Temperature  in  degrees  absolute  (Kelvin) 

TQ  *  Ice  point  temperature  (273.16°K)  at  standard  atmosphere 

T  *  Steam  point  temperature  (373.16°K)  at  standard 
5  atmosphere 

E  3  Saturation  vapor  pressure  of  pure  liquid  water  at 
s  steam  point  temperature  (1  standard  atmosphere  « 

1013.246  mb) 

EI0  »  Saturation  vapor  pressure  of  puie  water  ice  at  the  ice  point 
temperature  (0.0060273  standard  atmospheres  *  6.1077  mb) 

f  «  Correction  factor  for  the  departure  of  the  mixture  of  air  and 
water  vapor  from  the  ideal  gas  laws. 

*Ali  other  symbols  refer  to  pure  phase 
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I.  INTRODUCTION 

In  the  many  numerical  models  of  atmospheric  phenomena 
which  take  into  consideration  the  thermodynamics  of  moist  air, 
it  is  routinely  necessary  to  calculate  values  of  saturation 
vapor  pressure  one  or  more  times  at  each  grid  point  of  a  two 
or  three  dimensional  computational  domain  at  every  time  step. 
An  examination  of  computation  time  for  typical  models  would 
reveal  that  an  inordinate  amount  of  time  is  spent  in  the 
computation  of  vapor  pressure.  An  algorithm  for  computing 

this  parameter  with  accuracy,  and  which  requires  significantly 
less  computation  time  than  currently  used  procedures,  is 
presented  below. 


2.  SATURATION  VAPOR  PRESSURE  -  BACKGROUND 

As  Murray  ( 1 967 ) 1  pointed  out,  when  changes  in  phase  of 
water  substance  are  to  be  considered,  "...  the  saturation 
vapor  pressure  should  be  defined  as  the  equilibrium  vapor 
pressure  ..."  This  quantity  is  a  function  not  only  of  ambient 
temperature,  but  also  of  the  probable  impurities  in  water  or 
ice  substance,  and  the  size  and  form  of  droplets  and  crystals. 

In  practical  situations,  only  the  ambient  temperatures  can 
be  determined  with  any  degree  of  reality.  This  precludes  the 
use  of  equilibrium  vapor  pressure  in  thermodynamic  calculations. 
For  most  applications  (particularly  those  not  involving  explicit 
microphysical  calculation),  it  is  satisfactory  to  employ  satura¬ 
tion  vapor  pressure  (with  reference  to  plane  surfaces  of  water 
and  ice).  An  additional  complication  arises  from  the  difference 


Murray,  F.  W.,  1967:  On  the  computation  of  saturation 
vapor  pressure.  J .  Appl .  Meteor. ,  Vo 1 .  6,  pp.  203-204. 
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between  saturation  vapor  pressure  of  moist  air  and  the 
saturation  vapor  pressure  of  water  vapor  in  the  pure  phase. 
Saturation  vapor  pressure  of  the  pure  phase  is  a  function  of 
temperature  only.  Saturation  vapor  pressure  of  moist  air  is 
given  by 


(List,  1958)2 


This  can  be  contrasted  with  a  similar  formulation  for  the 
pure  phase 


(List,  1958) 


where  fw  is  the  ratio  of  the  saturation  vapor  pressure  for 
the  pure  phase  to  that  for  moist  air.  This  difference  arises 
from  three  related  effects  which  have  been  noted  by  List  (1958). 
This  difference  amounts  to,  at  most,  0.5%  in  the  ranges  of 
meteorological  pressure  and  temperature. 


3.  THE  STANDARD  -  GOFF-GRATCH  FORMULAS 

The  Goff-Gratch  formulations  are  the  internationally 
accepted  procedures  for  computing  saturation  vapor  pressure 
(in  the  pure  phase)  over  plane  surfaces  of  water  and  ice  (12th 
Conference  of  Directors  of  IMO,  Wash  ,  D.C.,  1947).  These 
formulations  are  as  follows: 


log10Ew  a  -  7.90298(TS/T  -1  )  +  5.02808  log1(J(Ts/T) 

-  1 ,3816X10~7(1011  *344(1-  T/Ts)_-|)  for  water  (3.1) 

+  8.1  328X10"3(10'3,491  49^Ts/T  -1)-1)  +  1  °g-|  qEws 


2 

List,  R.J.,  1958,  Smithsonian  Meteorological  Tables 
(6th  revision).  Washington,  D.C.,  The  Smithsonian  Institution, 
pp.  527. 
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log^j  3  -9.09718(TQ/T  -1  )  -  3.56654  log10(To/T) 


+  0.876793(1-  T/T0)  +  1og1()EI0>  for  ice  (3.2) 

where: 

Ew  =  saturation  vapor  pressure  over  a  plane  surface  of 
pure  ordinary  liquid  water  (mb), 

Ej  3  saturation  vapor  pressure  of  a  plane  surface  of 
pure  ordinary  water  ice  (mb), 

T  =  absolute  (thermodynamic)  temperature  (°K), 

T$  =  steam-point  temperature  (373.16°K), 

Ews  *  saturation  pressure  of  pure  ordinary  liquid  water 
at  steam-point  temperature  (  1  standard  atmos¬ 
phere  =  1013.246  mb) , 

EI0  =  saturation  pressure  of  pure  ordinary  water  Ice 
at  ice-point  temperature  (0.0060273  standard 
atmosphere  3  6.10/1  mb).  (Goff  and  Gratch,  1 946 )3 

The  forms  of  the  functions  above  are  quite  unwieldy  for  com¬ 
putational  purpose.  Murray  (1967)  transformed  them  through 
simple  operations  into  the  following. 


Ew  3  7.95357242X10 
w 
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(3.3) 


-  70242.1852  exp 
+  58.0691913  exp 


+  5.02808  In 
[  -8.03945282  M]|  ,  (3 


X  exp  j-18. 1  972839  |  J  +  5.02808  In 

r-26. 1205253 


3 

Goff,  J.A.  and  3.  Gratch,  1946:  Low-pressure  properties 
of  water  from  -110  to  212F.  Trans.  Amer.  Soc.  Heat.  Vent.  Enq.. 
fr^T-721. - 


E  j  =  5.75185606X10  exp  -20.947031 


(W 


3.56654  1 


nfM  2-01 
n \rr  tt 


01889049 

TOT] 


(3.4) 


Goff  and  Gratch  (1946)  claimed  a  2X10"^  percent 
uncertainty  for  the  water  reference  formulation,  Eq.  (3.1), 
(above  0°C)  and  3X10  percent  for  the  ice  reference  formula¬ 
tion,  Eq.  (3.2).  The  uncertainty  value  for  water  reference 
does  not  apply  to  values  below  0°C  where  no  experimental  data 
were  available.  Values  in  the  range  of  0  to  -50°C  were 
derived  by  direct  extrapolation. 

Murray's  reformulations,  Eqs,  (3.3)  and  (3.4),  of  the 
Goff-Gratch  formulas  differed  from  the  original  Goff-Gratch 
Eqs.  (3.1)  and  (3.2)  by  a  maximum  of  6X 1 0 “ 5  percent  (at  -25°C) 
and  3X1 0“ 5  percent  (at  -20°C)  for  water  and  ice  references 
respectively. 

4.  TETENS'  FORMULAS 

Murray's  transformations  of  the  original  Goff-Gratch 
equations  are  still  rather  unwieldy  from  the  standpoint  of  ease 
of  computation.  In  this  respect,  the  transforms,  Eqs.  (3.3) 
and  (3.4),  gain  little  over  the  original.  A  simpler  formulation 
for  determining  these  values  is  highly  desirable.  The  formula¬ 
tion  most  used  in  the  field  of  meteorology  has  been  and  is  that 
of  Magnus.  Tetdns  (1930)  gave  this  as 


i  F  s  M  +  w 

'10“s  t+v  w 


(4.1) 


^Tetens,  0.,  1930:  Uber  einige  meteorol ogi sche  Begriffe. 
)phys . ,  6,  297-309. 
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where  t  is  the  temperature  (°C),  w  =  0.7858  for  vapor  pressure 
in  mb,  and 

u  '  9,5  for  ice;  and  u  ~  7,5  over  water. 

v  =  265.5  v  =  237.3 

A  later  statement  of  this  formula  can  be  found  in  Haurwitz 
( 1 945 ) 5  and  is  given  as 

(— ) 

E$  *  6.1078  X  10  t+V  (4.2) 

where  u,  v  and  t  have  the  values  and  meanings  given  above. 

Murray  (1967),  for  the  purpose  of  achieving  greater  ease 
and  speed  of  computation,  reformulated  Eq.  (4.2)  to 


| 

55. 


I 

f 


i 

k' 


E$  =  6.1078  exp  (4.3) 

where  T  is  temperature  (°K),  and 

a-  21.8745584  f#r  ice.  a-1 7.2693882  for  water. 
b  =  7.66  b  =  35.86 

Murray  (1967)  showed  that  the  maximum  difference  between  the 
Goff-Gratch  and  Tetens  formulation,  for  both  ice  and  water, 
was  well  within  the  degree  of  uncertainty  demonstrated  by 
Goff  and  Gratch  (1946).  The  amount  of  error  (or  difference) 
arising  from  the  use  of  TetSns  formulation,  Eq.  (4.3),  can  be 
seen  in  Table  1  for  water  and  Table  2  for  ice.  The  maximum 
er"or  (difference)  for  water  is  4.4  percent  at  -50°C,  and, 
for  ice  is  3.0  percent  at  -50°C. 


'■'Haurwitz,  B.,  1945,  Dynamic  Meteorology,  New  York,  N.  Y., 
McGraw-Hill,  pp. 
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5.  TABLE  LOOK-UP  PROCEDURE 

In  addition  to  the  procedures  discussed  in  the  preceding 
sections,  there  is  another  method  for  computing  saturation 
vapor  pressure  which  is  quite  popular.  This  is  the  method  of 
"table  look-up."  It  is  particularly  attractive  when  consider¬ 
able  computer  memory  is  available  to  a  programmer.  The  "table 
look-up"  procedure  requires  the  storage  of  tabular  values  of 
saturation  vapor  pressure  over  a  desired  range  of  temperature. 
This  stored  table  becomes  a  permanent  and  integral  part  of  the 
program. 

For  a  given  temperature,  limits  of  saturation  vapor 
pressure  are  chosen  from  the  table.  This  is  accomplished  by 
determining  the  algebraically  largest  tabular  value  of  tempera¬ 
ture  less  than  the  temperature  in  question.  The  saturation 
vapor  pressures  for  this  and  the  next  higher  tabular  entries 
are  chosen  for  limits  (e.g.,  in  a  table  of  vapor  pressure 
values  for  each  whole  degree  temperature,  for  a  temperature  of 
6.55°C,  the  limits  will  be  the  vapor  pressure  values  f:  6°C 
and  7°C).  The  required  value  of  vapor  pressure  is  then  deter¬ 
mined  by  linear  (or  higher  order)  interpolation  within  these 
limits.  Linear  interpolation  is  normally  sufficient  because 
values  of  vapor  pressure  between  those  for  integer  values  of 
temperature  are  closely  approximated  by  a  straight  line. 

The  accuracy  attained  by  this  procedure  is  very  acceptable 
(see  Table  1)  with  the  largest  error  occurring  at  the  middle 
of  a  tabular  interval  (i.e.,  at  (T  +  0.5)°C  for  a  table  of 
values  at  1°C  tabular  intervals).  The  "table  look-up"  method 
is  more  accurate  than  the  Tetens'  formulations  for  temperatures 
less  than  -5°C  and  usually  less  accurate  for  those  above  -5°C 
(using  Goff-Gratch  as  the  standard).  It  is  slightly  more 
accurate  than  the  polynomial  procedures  (see  section  6)  from 
-50  to  -25°C  but  less  accurate  above  -25°C. 
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The  major  disadvantages  of  this  procedure  are  the 
requirement  for  considerable  memory  storage  and  the  necessity 
to  compute  interpolation  limits  which  is  time  consuming  (see 
Table  3). 


6.  A  POLYNOMIAL  APPROXIMATION 

Even  though  the  form  of  Tetens'  formula  given  by  Eq .  (4.3) 
and  the  procedure  of  table  look-up  are  less  complicated  than 
the  standard,  they  still  leave  something  to  be  desired  with 
respect  to  speed  of  computation.  This  is  so  because  of  the 
presence  of  the  exponential  function  in  Tetfe'ns  formulation 
which  is  quite  time  consuming  even  in  those  computers  having 
a  hard  wired  exponential  algorithm  (such  as  the  CDC-6500). 

The  table  look-up  procedure  also  consumes  considerable  time 
by  requiring  the  determination  of  interpolation  limits. 

Further,  the  procedure  requires  considerable  memory  storage. 

Clearly  an  even  speedier  procedure  for  the  determination  of 
saturation  vapor  pressure  is  highly  desirable. 

An  examination  of  the  curves  for  saturation  vapor  pressure 
indicated,  because  of  the  smooth  variability  with  changing 
temperature,  that  they  might  be  amenable  to  approximation  by 
some  order  of  polynomial  without  too  much  loss  of  accuracy. 

A  least  squares  fitting  procedure  was  applied  to  the  data 
for  the  temperature  range  from  -50°C  to  +50°C  for  the  water 
reference  and  -50°C  to  0°C  for  the  ice  reference.  Polynomials 
of  order  one  through  five  gave  a  very  poor  fit  with  maximum 
percentage  errors  of 


18212 

9294 

2854 

530 


for  order  1 


All  of  these  percentage  errors  occurred  at  -50°C  (for  the 
water  reference).  The  maximum  percentage  errors  for  the  ice 
reference  for  polynomials  of  order  one  through  five  are 

1540  for  order  1 

601  "  2 

149  "3 

26  "  4 

2.6  "  5  . 


The  sixth  order  polynomials  for  both  the  ice  and  liquid 
water  reference  gave  errors  of  less  than  one  percent  for 
the  entire  meteorological  range  of  interest.  The  polynomial 
formulation  for  saturation  vapor  pressure  is 

Es  3  a0  +  t  (a1+t(a2+t(a3+t(a4+t(a5+a6t))))),  (6.1) 

where  t  is  temperature  in  degrees  centigrade*  and  the  constants 
have  the  following  values 

for  water 

aQ  =  6.107799961  aQ  = 

a1  =  4 . 43651 8521 XI  O'"1  a]  = 

a2  =  1 . 428945805X1 0"2  a2  = 

a3  *  2 . 650648471 XI 0-4  a3  = 

a4  =  3.031 240396X1 0”6  a4  = 

a5  =  2 . 034080948X1 0“8  ag  = 

a6  »  6 . 1 36820929X1 0_1 1  ag  = 

(Range  of  validity:  -50°C  to  +50°C  for 
for  ice) 


for  ice 

6.109177956 
5 .034698970X1 0_1 
1 . 88601 3408X1 0"2 
4.176223716X10"4 
5.824720280X1 0"6 
4 . 8388031 74X1 0"8 
1 . 838826904X1 0"1 0 

water,  -5CC  +  0°C 


* 

The  coefficients  can  be  readily  re-evaluated  for  use 
with  temperatures  in  degrees  Kelvin. 
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7.  ACCURACY  RELATIVE  TO  THE  ACCEPTED  STANDARD 

Table  1  gives  values  of  saturation  vapor  pressure  (over 
water)  as  calculated  by  the  Goff-Gratch  formulation  (Eg), 

Eq.  (3.1);  by  Tetens'  formula  (Ey),  Eq.  (4.3);  by  table  look-up 
(ETL);  and  by  the  polynomial  (E^) ,  Eq.  (5.1).  Also  shown  in 
Table  1  are  values  of  percentage  departure  (error)  of  the 
Tetens  formula  and  polynomial  results  from  the  Goff-Gratch 
standard.  These  percentages  are  indicated  by  the  values  in 
parentheses.  Table  2  gives  analogous  information  for  satura¬ 
tion  vapor  pressure  with  respect  to  a  plane  ice  surface.  A 
quick  examination  of  these  tables  indicates  that,  with  only 
one  exception  (0°C,  for  ice,  Table  2),  the  percentage  departure 
due  to  the  polynomial  procedure  is  everywhere  many  times  less 
than  that  due  to  the  Tetens  formulation.  Therefore,  as  Murray 
(1967)  has  shown,  saturation  vapor  pressure  values  determined 
by  TetSns  formula,  Eq.  (4.3),  depart  from  the  standard  by  amounts 
less  than  the  degree  of  uncertainty  embodied  in  the  standard. 

The  polynomial  values  which  have  smaller  departures  must  be 
even  further  within  the  zone  of  uncertainty. 

Having  shown  that  the  polynomial  yields  values  of  satura¬ 
tion  vapor  pressure  which  are  at  least  as  accurate  as  the 
Tetens  formulations,  it  is  next  necessary  to  inquire  into  the 
relative  speeds  of  computation  of  the  methods  discussed  above. 

Each  of  the  procedures  was  used  to  compute  a  set  of  10,000 
saturation  vapor  pressures.  Evaluations  were  made  on  two 
computer  systems  --  the  CDC  3100  and  the  CDC  7600.  The  results 
of  these  evaluations  are  shown  in  Table  3. 

From  Table  3,  it  can  be  seen  that  the  polynomial  formula¬ 
tion  is  approximately  2.5  times  faster  than  the  best  Tetfe’ns 
formulation.  It  would  seem,  then,  that  the  demonstrated 
accuracy  and  speed  would  justify  the  use  of  the  polynomial  for 
the  determination  of  saturation  vapor  pressure  in  numerical 
model s . 
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Tabic  1.  Values  of  saturation  vapor  pressure  (in  mb)  (over 
water).  (Values  in  ()  are  ft  error.) 


E? 


7CC) 

£6 

£t 

*  £tl 

£l 

-50 

0.06354 

0.06078(-4.4) 

.0674(0.15) 

0. 06337 ( -0 . 29 ) 

-45 

0.11114 

0. 10735 ( -3.4) 

.1175(0.085) 

0.111 70 { +0 . 54 ) 

•40 

0.18914 

0.1 8423 ( -2.6) 

.1994(0.10) 

0.18915(+0.03) 

-35 

0.31387 

0.30783{-1 .9) 

.3301(0.12) 

0.31 3 1 4  ( -0 . 24 ) 

-30 

0, 50880 

0 . 501 77 ( -1 .38) 

.5338(0.10) 

0.50777(-0.20) 

-25 

0.80697 

0. 79928 (-0.95) 

.8448(0.09)  . 

0.80620(-0.10) 

•20 

1.25401 

1 ,24622(-0.62) 

1.3102(0.08) 

1 . 25386 (-0.01  ) 

-15 

1.91178 

1 . 90463 ( -0.37) 

1.9936(0.08) 

1 . 91 226 (+0. 02 ) 

-10 

2.86270 

2 . 85709 { -0.196) 

2.9799(0.07) 

2.86350(+0.028) 

-5 

4.21485 

4.21168(-0.075) 

4.3800(0.07) 

4.21548(+0.016) 

0 

6.1078 

6.1078(0.00) 

6.3370(0.06) 

6.1078(0.00) 

♦5 

8.71922 

8.72272(+0.04) 

9.0328(0.05) 

8.71839(-0.009) 

♦  10 

12.2723 

12.2  789 (♦0.054) 

12.696(0.04) 

1 2.2707 (-0.010) 

♦  15 

17.0438 

17.0S23(+0.050) 

17  608(0.05) 

17.0419(-0.012) 

♦20 

23.3730 

23. 381 0 (+0 . 034 ) 

24.117(0.04) 

23.3712 (-0.008) 

♦25 

31.6709 

31.6749(+0.013) 

32.640(0.04) 

31.6693(-0.005) 

♦  30 

42.4304 

42.4264(-0.009) 

43.678(0.04) 

42.4289(-0.002) 

♦35 

56.2366 

S6.2206(-0.028) 

57.829(0.03) 

56 . 2348 ( -0 . 002 ) 

♦40 

73.7775 

73. 7473 ( -0 . 041  ) 

75.790(0.03) 

73 .7754 ( -0 . 002 ) 

♦45 

95.8548 

95.8125(-0.044) 

95.8553(+0.0003) 

♦50 

123.400 

123.351 (-0.036) 

123.408(+0.006) 

•Temperature  for  this  column  are  offset  upwards  by  1/2  degree,  l.e.,  the  0JC 
entry  Is  really  the  value  for  ♦0.5°C. 


I’ 


* 

| 


14 


Table  2.  Values  of  saturation  vapor  pressure  (mb)  (over  ice) 
(Values  in  (  )  are  %  error.) 


T(°C) 

eg 

et 

el 

-50° 

0.03935 

0.03317 (-3.0) 

0.03963(+0.713) 

-45 

0.07198 

0.07032(-2.3) 

0 . 071 92 ( -0  -  076  ) 

-40 

0.12832 

0. 12611 (-1.7) 

0.12830(-0.002) 

-35 

0.22329 

0 . 22048 (-1.26) 

0. 2231 5 ( -0 . 066 ) 

-30 

0.37980 

0. 37644(-0 .89) 

0.37964 ( -0 . 041  ) 

-25 

0.63233 

0. 62857( -0 . 59) 

0.63230(-0.000) 

-20 

1.03173 

1 .02790(-0.37) 

1 .031 74 ( -0 .026) 

-15 

1.65167 

1 . 64824 ( -0 . 21 ) 

1.65141 (-0.035) 

-10 

2.59702 

2 . 59456 ( -0 . 095 ) 

2 .59646 (-0.021  ) 

-5 

4.01465 

4.01358(-0.024) 

4 .01462 ( -0 . 009 ) 

0 

6.10710 

6.1078(+0.011 ) 

6 . 1 091 8(  +  0. 035  ) 

Table  3.  Comparison  of  computation  times  for  SVP  procedures 
(over  water)  (units  -  sec/computation) . 


Procedure 

Eg  (3.3) 

ET  (4.2) 
'l 

Et  (4.3) 


3100 


16.316X10' 

9.849X10' 

4.583X10' 

4.975X10' 


7600 


20.331X10' 

11.400X10' 


5.76X10' 


(6.1) 


1.926X10' 


2.400X10' 


8.  COMPUTATION  OF  THE  DERIVATIVE  WITH  RESPECT  TO  TEMPERATURE 

Many  thermodynamic  computations  necessary  for  atmospheric 
simulation  require  determination  of  values  of  the  derivative 
of  saturation  vapor  pressure  with  respect  to  temperature. 
Differentiation  of  the  Goff-Gratch  equations  as  reformulated 
by  Murray  (see  Eqs.  3.3  and  3.4)  yields 


dEw 

W 


dEi 

TT 


-  j^5 . 02808  -  ^18. 


1973  +  446.844  exp  <-8.03945 


1834762  ....  -26.1205 

TTJTr  exp  {T7TT- 


!■ 

(  -26.1205  tlEw 

PW  F 


(8.1) 


for  liquid 


20.947(^0  +  3.56654  -  fj^yy 


T“ 


for  ice  (8.2) 


(Murray  &  Hollinden,  1966) 


It  can  be  seen  that  these  expressions  are  much  more  complicated 
than  even  the  original  Goff-Gratch  equations.  Besides  the 
complication  of  form,  it  is  also  required  to  calculate  the 
saturation  vapor  pressure  itself,  if  it  is  not  already  known  -- 
which  is  not  likely. 

Logarithmic  differentiation  of  Tetens'  formula,  Eq.  (4.3), 

gi  ves 

(8.3) 

=  5807.71  over  ice  and  4098.03 
over  water 

*  35.86  over  water  and  7.66  over 
ice . 


dEs 

ST 


A’Es 

rr%) 


2  where  A 


^Murray,  F.W.  and  A.B.  Hollinden,  1966:  The  evaluation 
of  cumulus  clouds:  A  numerical  simulation  and  'Its  comparison 
against  observations.  Douglas  Aircraft  Co.  Rep.  #SM-49372 . 
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This  expression  would  lead  to  a  very  rapid  calculation  of  the 
derivative  If  E$  Is  known.  If  It  Is  not,  then  the  calculation 
will  be  slightly  more  lengthy  than  Tetfe‘ns‘  calculation  for 
saturation  vapor  pressure. 

An  attempt  was  made  to  fit  a  polynomial  to  values  of  the 
derivative  over  the  range  of  temperatures  of  meteorological 
Interest.  The  data  used  was  obtained  by  evaluations  of  Eqs. 
(8.1)  and  (8.2).  Polynomials  of  order  5,  6  and  7  all  showed 
acceptable  error  patterns  (l.e.,  errors  less  than  those  arising 
from  the  TetSns  formulation,  Eq.  (8.3))  for  water.  For  Ice 
reference,  orders  5  and  6  were  acceptable,  but  strangely  enough 
higher  orders  were  not.  An  evaluation  of  mean,  maximum  and 
root-mean-square  errors  indicated  that  the  sixth  order  polyno¬ 
mial  was  again  the  optimum  choice.  The  polynomials  take  the 
same  form  as  Eq.  (6.1)  but  with  coefficients  as  shown  below: 


for  water 

aQ  -  4.438099984X10"1 
a1  •  2.857002636X10"2 
a ^  »  7. 938054040X10"4 
a3  -  1 . 21 521 5065X1 0”5 
a4  *  1 . 036561403XT0”7 
a5  *  3. 53242181 0X1 0"1 0 
a6  *  -7 . 090244804X1 0"1  3 

(range  of  validity: 


for  Ice 

aQ  =  5.030306237X10"1 
a1  »  3.773255020X1  O’2 
a2  «  1 .267995369X10’3 
a3  =  2 . 4775631 08X1 0"5 
a4  =  3 . 0056931 32X1 0”7 
a5  =  2 . 1 58542548X1 0”9 
a6  =  7.131 097725X1 0 _1 2 

-50°C  to  +50°C  for  water; 

-50°C  to  0°£  for  i ce) 


(8.4) 
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Table  4  shows  the  values  of  ^4  as  computed  using  the  Goff- 
6ratch  { Dq )  derivatives,  TetSns1  derivative  (DT),  and  those 
determined  from  the  polynomial  (DL)  just  discussed.  Table  5 
shows  similar  Information  for  dE  /dT  for  the  ice  reference 
case.  As  with  the  primary  functions,  the  polynomial  expres¬ 
sions,  coefficient  set  (8.4),  for  the  derivatives  (with 
reference  to  both  ice  and  water)  show  departures  (errors) 
from  the  Goff-Gratch  standard  which  are  considerably  less  than 
those  determined  by  the  use  of  TetSns'  formulation,  Eq.  (8.3). 
Derivative  computation  times  are  comparable  (as  would  be 
expected)  to  those  for  the  primary  expressions  for  saturation 
vapor  pressure  (see  Table  3). 

9.  SUMMARY 

It  has  been  shown  that  it  is  possible  to  formulate  a 
polynomial  approximation  for  both  saturation  vapor  pressure 
and  its  derivative  with  respect  to  temperature  that  is  at 
least  and  's  generally  much  more  accurate  than  currently  used 
procedures  (Tet6ns'  formula).  Accuracy  was  measured  in  terms 
of  departure  from  values  derived  from  the  Goff-Gratch  formulas 
which  are  the  internationally  accepted  standards.  The  poly¬ 
nomial  errors  are  well  within  the  degree  of  uncertainty 
connected  with  the  Goff-Gratch  procedures.  The  polynomial 
procedures  have  been  demonstrated  to  consume  significantly 
less  computer  time  than  methods  currently  in  use.  The  employ¬ 
ment  of  this  procedure  will  result  in  significant  savings  in 
the  consumption  of  computer  resources  and  money.  These  same 
polynomials  may  be  used  to  evaluate  actual  vapor  pressure  by 
using  the  dew  point  temperature  in  lieu  of  air  temperature. 


Table  4.  Derivative  of 
(mb/degree) . 


saturation  vapor  pressure  over  water 


T{°C) 

°G 

r  ®t 

Di 

-50°C 

0.007286 

0.007100(-2. 6) 

0.007188(-1 . 35) 

-45°C 

0.012113 

0.011 897 ( -1 .8) 

0 .01 2234 (+1 . 001 ) 

-40°C 

0.019624 

0.01 9394 ( -1 .17) 

0.019644(+0.099) 

-35°C 

0.031042 

0.030824(-0.70) 

0 . 030940 ( -0 . 329) 

-30°C 

0.048021 

0 . 047849 (-0.36) 

0.047887 ( -0 .279) 

-25°C 

0.072756 

0.072673(-0.11 ) 

0 . 072678 ( -0.107) 

-20°C 

0.10811 

0. 10816(+0.045) 

0 . 1 0812 (+0.013) 

-1 5°C 

0.15773 

0. 15795(+0.1 36) 

0 . 1 5782 (+0 . 059 ) 

-10°C 

0.22622 

0. 22662 (+0 .176) 

0 . 22634 (  + . 052 ) 

-5°C 

0.31927 

0.31 983(+0 .179) 

0 . 31 935 (+ . 025 ) 

0°C 

0.44381 

0.44449(+ . 1 54) 

0.44381 (0.00) 

5°C 

0.60817 

0.60886(0.114) 

0 . 60809 (-0.013) 

10°C 

0.82225 

0.82279(0.065) 

0.82211 (-0.016) 

1 5°C 

1.0976 

1.0978(0.016) 

1 . 0975 ( -0.010) 

20°C 

1.4477 

1 . 4473 ( -0 . 028) 

1 . 4476 ( -0 . 005 ) 

24°C 

1.8878 

1  -  8867 ( -0 . 063 ) 

1 .8878(+0.002) 

30°C 

2.4354 

2.4334(-0.083) 

2 . 4355 (  +  0 . 003 ) 

35°C 

3.1100 

3 . 1 072 ( -0 . 087 ) 

3 . 1 1 00 ( -0 . 0002 ) 

40°C 

3.9331 

3, 9303 ( -0 . 072 ) 

3.9331 (-.0009) 

45°C 

4.9287 

4. 9269  ( -0 . 036 ) 

4 . 9286 ( -0 . 001 ) 

50°C 

6.1228 

6.1241 (+.022) 

6 . 1 230 (+ . 003) 
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Table  5.  Derivative  of  saturation  vapor  pressure  over  ice 
(mb/degree) . 


T('C) 

°G 

°T 

Di 

-50°C 

0.004859 

0 . 004773 ( -1 .8) 

0.004875(+0.32) 

-45°C 

0.008506 

0.008400(-l .23) 

0.008502 ( -0 . 033) 

-40°C 

0.014521 

0.014403(-0.81 ) 

0.014512(-0.059) 

-35°C 

0.024219 

0.024101 (-0.49) 

0.024210(-0.036) 

-30°C 

0.039519 

0.039420(-0.25) 

0.039515 (-0.010) 

-25°C 

0.063166 

0.0631 1 5 ( -0 . 081  ) 

0.063166(-0.000) 

-20°C 

0.099022 

0.099050(+0.028) 

0.09901 3( -0. 009) 

-1  5°C 

0.15241 

0.1 5255 (+0 . 089 ) 

0.152381 (-0.018) 

-10°C 

0.23058 

0.23083(+0.109) 

0.23053(-0.023) 

0.34317 

0.34350(+0.097) 

0. 3431 5 ( -0 . 005 ) 

0.50292 

0. 50322 (+0. 061 ) 

0 . 50303 (+0 . 022 ) 
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APPENDIX 


SATURATION  VAPOR  PRESSURE  AT  VERY  LOW  TEMPERATURES 

The  polynomial  approximations  discussed  and  evaluated  in 
the  foregoing  were  strictly  valid  only  for  given  ranges  of 
temperature  (-50°C  to  +50°C  for  water  reference  and  -50°C  to 
0°C  for  ice  reference).  Attempts  to  use  these  polynomials 
at  lower  temperatures  resulted  in  drastic  departures  from  the 
Goff-Gratch  values.  To  overcome  this  difficulty,  values  of 
saturation  vapor  pressure  (with  respect  to  both  ice  and  water) 
were  generated  from  the  Goff-Gratch  formulas.  These  data  were 
then  used  to  derive  a  new  set  of  polynomial  approximations. 

The  polynomials  are  again  of  the  sixth  order  and  take  the  same 
form  as  that  in  Eq.  (6.1)  with  the  coefficients  shown  below: 


for  water 

a„  =  4.866786841 
o 

a]  =  3. 1 52625546X1 0"1 
a2  =  8.6401 88586X1 0”3 
a3  =  1 . 279669658X1 0"4 
a4  =  1 . 07795591 4X1 0"6 
a5  *  4. 886796102X1 0~9 
ag  =  9.296950850X10"  12 

(Range  of  validity: 


for  ice 

aQ  =  3.927659727 
a]  =  2 . 643578680X1 0"1 
a2  =  7 . 505070860X1 C~3 
a3  =  1 . 1 47668232X1 0~4 
a4  =  9 . 948650743X1 0"7 
a5  =  4. 626362556X1 0"9 
a6  =  9.001382935X1 0"1 2 

-100°C  to  ~50°C ) 


Table  A-l  shows  the  accuracy  of  the  approximation  with  respect 
to  water  as  compared  to  that  of  Tetfe’ns  using  Goff-Gratch  as 
a  standard.  Percentage  errors  are  shown  in  parentheses.  With 
the  exception  of  -100°C,  the  polynomial  values  depart  from  the 
Goff-Gratch  values  by  a  considerably  smaller  percentage  than 


Preceding  page  blank 


23 


those  derived  from  Tetens  formulation.  Table  A-2,  shows 
similar  information  for  saturation  vapor  pressure  with  respect 
to  ice  reference. 


Table  A-l.  Values  of  Saturation  Vapor  Pressure  (mb) 
(water  reference)  at  very  low  temperatures. 


t(°c) 

eg 

et 

el 

-100 

. 24X1 0"4 

.21X1 0~4 (-12 .25) 

.33X10"4(36.14) 

-95 

.70X10"4 

.  60X1 0_4( -1 4 .12) 

.68X1 0~4 (-2.49) 

-90 

. 1 87X1 0"3 

.160X10”3{ -14.41 ) 

.  1 92X1 0“3 (  +  3 . 1 3 ) 

-85 

. 462X1 0"3 

. 398X1 0  ~  3 ( - 1 3.77) 

. 458X1 0”3 ( -0.72) 

O 

CO 

3 

.  1 07X1 0-2 

.  936X1 0"3( -1 2 .60) 

. 1 06X1 0~2 ( -0.74) 

-75 

.235X1  O'2 

.  209X1 0"2 ( -1 1.17) 

. 235X1 0~2 (+0 . 00 ) 

-70 

. 492X/0"2 

. 444X1 0“2 (-9.65) 

. 492X1 0~2 (  +  0 . 07 ) 

-65 

. 985X1 0"2 

. 905X1 0"2 (-8.16) 

.984X1 O-2 ( -0 . 1 3 ) 

-60 

.01898 

.01 769( -6. 76) 

.01  894 ( -0.18) 

-55 

.03529 

.03335 { -5  .49 ) 

.03528(-0.03) 

-50 

0.06356 

. 06078 ( -4 . 37 ) 

. 06362 (+0.11 ) 

Table  A-2.  Values  of  Saturation  Vapor  Pressure  (mb)  (ice 
reference)  at  very  low  temperatures. 


t(°c) 

eg 

Er 

el 

O 

o 

1 

.  1 403X1 0'4 

.1114X1 0“4(-Z0. 61 ) 

.231 3X1 0"4 (+64 . 85 ) 

-95 

.  3784X1 0  ~4 

. 31 1 7X1 0-4 ( -1 7 . 62 ) 

.  3587X1 0"4 (-5.20) 

-90 

.  9672X1  0**4 

. 3224X1 0"4 (-14.96) 

.  1 047X1 0~3 (+8 . 23 ) 

-85 

.  235X1 0”3 

.  2056X1 0"3 ( -12.61  ) 

.2321 XI 0'3 ( -1 .39) 

o 

cc 

1 

.  5472X1 0"3 

.  4894X1 0~3( -10. 55) 

.  5375X1  0"3 ( -1  .78) 

-75 

.  1 220X1  O’"2 

.  1 i 1 3X1 0" 2( -8 . 75 ) 

.  1 222X1 0’2 (  +  0 . 02 ) 

-70 

.  261 5X1 0”2 

.  2427X1 0’2 ( -7 . 1 9 ) 

.  2622X1 0"2 (  +  0.29) 

-65 

.  5406X1 0'2 
.  1 080X1 0”1 
.2092X1 0'1 
.  3935X1 0"1 

.  5090X1 0"2 (-5.84) 

.  1030X1 0"1 (-4.69) 

.  201 5X1 0"1 (-3.71 ) 

.  5396X1 0"2 (-0.17) 

-60 

.  1077X10'1 (-0.32) 

-55 

.  2091  XI 0"1 (-0.06) 

-50 

.  3821  XI 0"1  (-2.89) 

.  3940X1 0”1 (  +  0.15) 
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The  polynomial  approximations  for  the  derivative  of 
saturation  vapor  pressure  with  respect  to  temperature  for  ice 
and  water  reference  is  also  of  sixth  order  and  have  the  same 
form  as  Eq.  (6.1). 


for  water 

a„  -  4.086240791X10 
a1  =  2.51611 8369X10 
a2  =  6.576862688X10' 
a3  =  9 . 325531 513X1  O' 
a„  =  7.550718726X10 


for  ice 


-1 


3.303373957X10 

6.088242842X10 


-8 

-10 

•13 


-1 

-2 

-4 


a„  =  3.808584154X10 
a1  =  2.482536693X10 
a2  =  6.847888228X10 
a3  =  1.020750000X10 
a4  =  8. 651 642035X1 0~8 


-5 


a5  =  3.944818571X10 
ag  =  7.544801596X10 


-10 

-13 


(Range  of  validity:  -100°C  to  -50°C) 


Tables  A-3  and  A-4  show  the  comparative  accuracy  of 
polynomial  approximations  and  Tetfe'ns'  formula  for  the  deriva¬ 
tive  of  vapor  pressure  with  respect  to  water  and  ice  reference. 
It  is  interesting  to  compare  the  juncture  (-50°C)  of  the  tables 
in  this  appendix  with  the  corresponding  tables  in  the  main 
body  of  the  text  (Table  A-l  corresponds  to  Table  1;  A-2  to  2; 
A-3  to  4;  and,  A-4  to  5).  For  vapor  pressure  with  respect  to 
water,  the  "normal"  value  (normal  implies  an  ordinary  meteor¬ 
ological  temperature  range  as  opposed  to  the  very  low 
temperature  range)  is  .06337  mb  (-.29)  while  the  low  range 
value  is  .06362  mb  (+0.11).  The  absolute  percentage  difference 
between  these  values  is  0.39%  with  the  low  range  value  being 
the  more  accurate  with  respect  to  the  Goff-Gratch  values.  For 
vapor  pressure  with  respect  to  ice,  this  percentage  is  2.58%; 
for  the  derivative  with  reference  to  water,  1.45%;  and  lastly, 
for  the  derivative  with  reference  to  ice,  0.23%.  In  all  cases. 


Table  A-3.  Values  of  derivative  of  saturation  vapor  pressure 
(mb/deg)  with  respect  to  temperature  (water  reference) 
at  very  low  temperatures. 


t(°c) 

dg 

dt 

ci 

-100 

.  538X1 0"5 

.458X1  0”5(-l 4.91 ) 

.  593X1 0"5 (  +  1 0 .28) 

-95 

.  1 432X1 0"4 

.121 6X10-4(-l 5. 05) 

.1421 XI 0"4(-0 . 75 ) 

-90 

.  351 4X1 0”4 

. 301 7X1 0"4 (-14.13) 

.3530X1 0"4 (+0 .47 ) 

-85 

.  8055X1 0"4 

.  7035X1 0"4 (-12.66) 

.  8029X1 0”4 (-0.32) 

-80 

.  1  742X1 0“3 

.  1  551 XI 0~3 ( -1 0 . 95 ) 

.  1 738X1 0~3 ( -0.22) 

-75 

.  3581 XI 0"3 

.  3251  XI 0"3( -9 . 20) 

. 3581 XI 0~3(+0 . 00 ) 

t 

O 

.  7037X1 0"3 

.  6507X1 0"3( -7 . 53) 

.  7037X1 0~3(-0 .00) 

-65 

.  1 329X1 0“2 

.  1 249X1 0“2 ( -6.00) 

. 1 327X1 0-2 (-0.09) 

-60 

.2419X10'2 

.  2307X1 0"2 (-4.65) 

.241 7X1 0"2 (-0.10) 

-55 

. 4262X1 0"2 

.  41 1 2X1 0“2 (-3.51) 

. 4262X1 0”2 (-0.02) 

-50 

. 7286X1 0”2 

.  71 00X1 0~2 ( -2.55) 

.  7291  XI 0"2 (  +  0 . 07 ) 

Table  A-4.  Values  of  derivative  of  saturation  vapor  pressure 
(mb/deg)  with  respect  to  temperature  (ice  reference)  at 
very  low  temperatures. 


t(°c) 

°G 

°T 

dl 

-100 

,  286X1 0'5 

.  236X1 0-5 ( -1  7.50) 

. 358X1 0-5 (+25 . 18) 

-95 

.  730X1 0-5 

.  623X1 0"5 ( -14.70) 

.  71 5X1 0~5(-2 .03) 

1 

lO 

o 

.  1 766X1 0"4 

.  1 550X1 0”4 ( -1 2.23) 

. 1 81 9X1 0-4 (+2 . 95 ) 

1 

CO 

cn 

. 4076X1 0”4 

.  3665X1 0"4 ( -10.07) 

.4048X10”4(-0.67) 

-80 

. 8997X1 0"4 

. 8259X1 0"4 ( -8.20) 

.  8929X1 0-4 (-0.76) 

-75 

.  1 907X1 0-3 

.  1 781 XI 0"3( -6 . 60 ) 

.  1907X10~3(  +  0.00) 

-70 

.  3891X1  O'3 

.  3687X1 0"3( -5.24) 

. 3895X1 0"3(+0. 09) 

-65 

.  7665X1 0-3 

.7352X10"3(-4.08) 

.  7655X1 0  3 ( -0.1  3) 

-60 

.1462X1 0-2 

.  1 41 6X1 0"2 ( -3 . 1 2 ) 

.  1 459X1 0”2 (  -0 . 1 9 ) 

-55 

.2703X1  O'2 

.  2640X1 0”2 ( -2 . 33 ) 

.  2702X1 0"2 ( -0 . 03 ) 

-50 

.  4859X1 0-2 

.  4777X1 0"2 ( -1  .68) 

.  4864X1 0 "2 (  +  0 . 1 1 ) 
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the  cold  range  value  is  the  more  accurate.  The  important 
point  to  make,  however,  Is  that  the  polynomials  for  both  the 
"normal"  and  the  very  low  ranges  may  be  used  together  without 
suffering  a  severe  or  significant  discontinuity  in  variation 
in  the. curve  of  values  (both  zero  and  first  order). 

Finally,  computation  time  for  the  various  procedures  are 
comparable  tc  those  shown  in  Table  3  and  maintain  the  same 
relativity. 


